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A newly identified Luminous Blue Variable in the galactic 
starburst cluster Westerlund 1 * 
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Abstract. We present observations of the massive transitional star W243 in Westerlund 1. We find an apparent 
spectral type of early-mid A from our data, in contrast to an earlier classification of B2I, made from data obtained 
in 1981. The concurrent development of a rich emission line spectrum suggests a very high mass loss rate; the 
continued presence of He I emission suggesting that the underlying star remains significantly hotter than implied 
by its apparent spectral type. We suggest that W243 is a Luminous Blue Variable undergoing an eruptive phase, 
leading to an extreme mass loss rate and the formation of a pseudo photosphere. 
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1. Introduction 

Before becoming hydrogen depleted Wolf Rayets (WR), 
massive stars pass through a period of enhanced mass loss, 
which removes the bulk of the H-rich mantle over their 
cores. However, passage through the 'zoo' of transitional 
objects - Red Supergiants (RSGs), Yellow Hypergiants 
(YHGs), Luminous Blue Variables (LBVs) and sgB[e] 
stars - is at present poorly understood, largely due to the 
rarity of such objects. Nevertheless, all such objects ap- 
pear to share certain characteristics, notably mass loss 
rates of >lO -5 M0yr _1 , significantly slower, denser winds 
than either their O star progenitors or WR progeny and 
highly variable surface temperatures and radii. Indeed, 
significant variability across all observable physical param- 
eters appears to be the defining characteristic of this brief 
episode of stellar evolution. Clearly, observations of such 
stars in clusters and their subsequent placement in HR 
diagrams - from which ages and progenitor masses could 
be inferred - would greatly enhance our understanding. 
Unfortunately, examples of such clusters are understand- 
ably rare and presently provide rather poor evolutionary 
constraints. 

Westerlund 1 (henceforth Wd 1; Westerlund I1961|) is 
a highly reddened cluster found at a distance of ~2.5kpc 
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(Clark et al. 2004). Prior to 2001, the only spectroscopic 
survey of Wd 1 had been that of Westerlund ( 19$?| hence- 
forth We87) which encompassed ~20 of the brightest clus- 
ter members. Motivated by an unusually rich population 
of radio sources within Wd 1 (Clark et al. 1199 8 Dougherty 
et al., in prep.) the cluster was subsequently spectroscop- 
ically re-observed in 2001, with follow up observations 
made in 2002-3. The resultant data revealed a hitherto un- 
suspected population of massive post-MS objects in Wd 1 
(Clark & Negueruela 123321 Negueruela & Clark 123351 and 
Clark et al. l2004f) . With a mass comparable to other star- 
burst clusters in the Local Group (e.g. Negueruela & Clark 
2003| and an estimated age of 3-5 Myrs, Wd 1 hosts a 
unique population of massive transitional objects. Here we 
present observations of one such star - W243 (RA=16 47 
07.5 <5=-45 52 28.5; J2000) - that demonstrate significant 
spectral variability, both internally and in comparison to 
the data of We87, suggesting a likely identification as an 
LBV undergoing a significant eruptive event. 



2. Observations & data reduction 

Observations of W243 have been made at a number of 
different facilities over the past two years; these are sum- 
marised in Table 1. Data reduction employed both pack- 
ages within the Starlink software suite and also MIDAS; 
further details are provided in Clark & Negueruela ( 2002 ) 
and Clark et al. (|2004[1 . Data obtained with the NTT and a 
subset of the VLT dataset are presented in Figs. 1-3; pre- 
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Date Instrumental Configuration(s) 
23/6/2001 ESO 1.52m, Loral #38 camera 

GRAT 1 (6000-llOOOA at ~5A/pixel) 

EW Htt = -18±4A 
7/6/2002 NTT + EMMI in REMD mode (2x2 binning) 

GRAT 7 (6310-7835A at ~0.8A/pixel) 

EW Htt = -18.7±0.3A 

GRAT 6 (8225-8900A at ~0.4A/pixel) 
6/6/2003 NTT+EMMI in REMD mode (2x2 binning) 

GRAT 6 (6440-7140A at ~0.4A/pixel) 

EW Htt = -20.7±0.3A 
7/6/2003 NTT + EMMI in REMD mode (no binning) 

GRAT 6 (8225-8900A at ~0.2A/pixel) 
8/6/2003 NTT + EMMI in RILD mode (no binning) 

GRISM 4 (6000-10700A at ~0.4A/pixel) 
21/9/2003 VLT UT2 + UVES (red arm) 

CD#4 (6600-10600A; R « 40 000 ) 

CD#3 (4760-6840A R w 40 000) 
EW Ha = 22.3+0.2A 

Table 1. Observation log for W243 between 2001-3, giv- 
ing telescope, instrument, instrumental configuration and 
resultant wavelength range and dispersion. Note the obser- 
vations on 8/6/03 were made through clouds. Note for rea- 
sons of brevity we have also listed the Equivalent Widths 
of the Ha profiles from each observation; due to the low 
resolution of the 2001 data (not shown) the error associ- 
ated with this measurement is substantially higher than 
subsequent observations. 
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Fig. 1. Left panel: Ha profile obtained with UVES 
(21/9/2003) - note the narrow (FWHM=114kms^ 1 ) dou- 
ble peaked profile (separation=48kms _1 , central absorp- 
tion blucshiftcd 31kms _1 from rest wavelength) and 
emission wings extending to projected velocities of 
±10 3 kms _1 . Right panel: Comparison of blended profiles 
of Pall and Ni 8747A from the NTT (2003 June 7; red 
dotted lines) and the VLT (2003 September 21; solid black 
lines, resolution degraded to match the NTT spectrum) 
demonstrating continued spectral evolution over only 3 
months. 



sentation and quantitative spectroscopic analysis of the 
full VLT dataset is defered for a future publication. 

3. Discussion 

The observations between 2001-3 reveal that W243 pos- 
sesses a complex spectrum with emission from H I (Ha & [3 
- Fig 1, latter not shown) and He I (6678A & 7065 A- Fig 2). 
The Ha profile has significantly strengthened during this 
period as indicated in Table 1. Weak Can emission is also 
present in 2002; by 2003 this emission had considerably 
strengthened (Fig. 3) while emission was also observed in 
O I 8448A. These chan ges arc also reflected in the stronger 
Paschen lines (Figs. 1 & 3), with an increase in the depth 
of the absorption feature between 2002-3 accompanied by 
an increase in strength of the emission component of the 
profile. Fig. 1 clearly indicates that these changes have 
continued in the 3 month period between the NTT and 
VLT observations; the Call and Oi 8448 A lines have like- 
wise increased in strength (not shown). 

Higher Paschen series lines are observed to be in ab- 
sorption as are a number of N I transitions longwards of 
~8560A. The O I 7774A complex is also found to be in 
absorption, while the Oi 8448A and Can lines demon- 
strate apparent inverse P Cygni profiles. While contam- 
ination of the absorption lines by wind emission compli- 
cates spectral classification - most notably for the Paschen 
series and Ca n - the presence of the strong N I permits a 
broad classification. Clark ct al. ( 2004 ) find that such lines 
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Fig. 2. Narrow (FWHM~ 60 ± 2kms~ 1 ), single peaked 
profiles of He I 6678A (left panel; EW=-0.24±0.04A) and 
7065A (right panel; EW=-0.78±0.04A) from UVES data 
(rest wavelengths indicated). 



are absent for supergiants with spectral types earlier than 
^B3-5, very weakly in absorption for B5-8 stars, show a 
rapid increase in strength between B8-A2 before remain- 
ing constant for later A stars. Comparison of our spectra 
to those of Munari & Tomasella l|1999l Fig. 3) show that 
W243 is clearly later than B8 and is consistent with an A2 
supergiant. Unfortunately, the strength of the Can lines is 
the primary discriminant for the A subtypes so we cannot 
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Fig. 3. I band (8200-9000 A) spectra of W243 obtained between 2002-2003 compared to spectra of B & A supergiant 
standards (data from Munari & Tomasella 1999). Note the overlap in wavelength ranges betwen upper and lower 
panels. Principal transitions are indicated; the broad absorption feature at ^8620 A is a DIB and is common to all 
our spectra of Wd 1 members. The spectrum of the peculiar YHG IRC+10 420 is also presented; a wavelength shift 
has been applied to this spectrum to correct for the systemic and intrinsic velocity shifts identified for this object by 
Oudmaijer H1998|l . 



exclude a later subtype; the lack of Fe I absorption lines - 
present in F supcrgiants - precluding a classification later 
than A. Note however that the presence of He I emission 
is clearly discrepant; we return to this below. 

Based on observations obtained in 1981, We87 re- 
ported a spectral type of B2I for W243 and a spectrum de- 
void of any emission features. Given the resolution of the 
spectra, this classification was based on the presence of 
weak Paschen absorption lines and an absence of Call ab- 
sorption features (Bengt Westerlund, priv. comm. 2003). 
As demonstrated by our spectra, wind contamination can 
render the Paschen and Ca n lines unsuitable for classifica- 
tion, leading to an erroneously early classification for low 
S/N and resolution data. However, consideration of real 
and synthetic spectra (e.g. Clark et al. 200|J) shows that 
the O I 7774A feature seen weakly in absorption in the 
1981 spectrum is consistent with an early-mid Bla spec- 



tral type. Since significant infilling of this feature is not 
observed in our spectra - despite the presence of numerous 
strong wind emission lines - we conclude that the weak- 
ness of this feature in the 1981 spectrum was not due to 
wind contamination, given that no strong wind lines (e.g. 
Ha) were present at that time. Hence, we conclude that 
differences in the absorption features in We87's and our 
spectra are not solely the result of a combination of wind 
contamination and differing spectral resolutions and that 
the classification of W243 as an early B supergiant in 1981 
is supportable. Our confidence is further strengthened by 
the excellent correspondence between the spectral classifi- 
cations of the remaining 6 B & A supergiants both studies 
have in common. 

Consequently, we are driven to the conclusion that be- 
tween 1981-2001 the apparent spectral type of W243 has 
changed from ~B2-5Ia to A2I (or later - implying a de- 



4 



J. S. Clark and I. Negueruela: A new LBV in Wd 1 



crease in temperature of >10kK); a change accompanied 
by the development of a rich emission line spectrum. 

Unfortunately we have no contemporaneous photome- 
try of W243, but adopting the values presented by We87 
and assuming that the evolution in spectral type occurred 
at constant bolometric luminosity, we may infer a lumi- 
nosity of ~6.3xl0 5 L© for W243 (distance and reddening 
from Clark et al. !20,(3J- At such a luminosity, evolution 
from spectral type B2I to A2I (or later) results in the 
star crossing the Humphreys-Davidson limit and passage 
into the so called 'yellow void' (e.g. de Jager[2001). Of 
particular interest therefore, is the remarkable similarity 
between W243 and the peculiar YHG IRC+10 420, which 
is thought to be rapidly evolving bluewards across the yel- 
low void (e.g. Humphreys et al. 1201)21 Oudmaijer [1998). 

A comparison of spectra is shown in Fig. 3 - to the 
best of our knowledge no other cool, luminous star demon- 
strates a comparable emission line spectrum 1 . We find 
striking similarities between the line profiles of the strong 
wind lines of H I & Ca n and the broad Ha emission wings 
of both stars. The emission wings likely result from elec- 
tron scattering and while the origin of the double peaked 
profiles is uncertain, recent observations of IRC+10 420 by 
Humphreys et al. ( 2002;) apparently reveal no large scale 
departure from spherical symmetry for the stellar wind 
(note that no other LBV is known to demonstrate a dou- 
ble peaked Ha line profile). Differences do exist between 
the 2 objects - notably in the presence of Fen emission 
in IRC+10 420 and the presence of strong O I 8446A and 
- in particular - Hei emission in W243. Moreover, we do 
not find the emission lines to be blue shifted with respect 
to the systemic velocity as Oudmaijer (1998) found for 
IRC+10 420. 

A further difference between the two stars is the pres- 
ence of a complex dusty ejection nebula surrounding IRC 
+10 420; which Humphreys et al. (|1997|) find to extend 
to ~ 5" (0.13 pc at a distance of 5kpc). Consideration 
of both a high resolution Ha image - obtained under ex- 
cellent (~ 0.4") seeing in 2003 June - and the Midcourse 
Space Experiment 8-25/im fields for W243 reveal no ev- 
idence for comparable ejecta. We consider two possibil- 
ities for this observation. Firstly, the nebula associated 
with IRC +10 420 is thought to have formed in a proceed- 
ing RSG phase. However, W243 has evolved from higher 
temperatures in the past 20 years, suggesting that such 
a nebula may yet form. Alternatively, the harsh environ- 
ment of Wd 1 is clearly inimical to the long term survival 
of dusty ejecta; the radiation fields and winds of several 
hundred massive cluster members may have destroyed any 
such nebula (although compact, presumably young, neb- 
ulae are associated with W9 & 26; Clark et al.[1998). 

Therefore, while W243 and IRC +10 420 may be in 
different evolutionary stages either side of an intervening 
RSG phase, we consider it likely that they both share sim- 

1 Albeit with differences in line strengths; Ha is significantly 
stronger in W243 than IRC +10 420 while the opposite is true 
for Call emission. 



ilar wind properties. For IRC +10 420, Humphreys et al. 
(2002) suggest that a slow wind with a particularly high 
mass loss rate (3-6 xlO - M©yr _1 ) is sufficiently dense to 
permit the formation of a cool (^8000K) 'pseudo' pho- 
tosphere which veils the underlying star (e.g. Davidson 
1987). Such an explanation for W243 would then naturally 
explain the apparent late spectral type, narrow emission 
lines, broad Ha emission wings and most importantly, the 
presence of He I emission, which implies an earlier spectral 
type than ^mid A. 

Dougherty et al. (in prep.) find time averaged ra- 
dio fluxes of 1.04±0.07mJy (8.64GHz) and 0.95+0. lOmJy 
(4.8GHz) for W243 between 2000-2. The spectral index is 
therefore flatter than expected for a canonical stellar wind, 
although still consistent with thermal emission. Adopting 
the wind properties given for the B8Ia star HD 160529 
(Leitherer et al. 1 995 note that the winds of A and later 
supergiants are poorly understood) for W243 yields a mass 
loss rate of 4.1xl0 -6 MQyr~ 1 . While unexpectedly small 
in light of the preceding discussion, it should be remem- 
bered that this estimate assumes that W243 has the stellar 
properties of a late B supergiant and that the wind is fully 
ionised. Clearly both assumptions may be in error; hence 
we choose to regard this value as merely a lower limit to 
the true mass loss rate - if the radio emission arises from 
the wind. 

We therefore suggest that W243 be considered as a 
new addition to the known galactic LBVs. Furthermore, 
we suggest that its present 'composite' spectral appear- 
ance is due to a continuing increase in mass loss rate in 
a LBV eruption, resulting in the formation of a 'pseudo' 
photosphere. Support for such a conclusion is provided by 
M33 Var B, which during such an eruption also displayed a 
late spectral type with anomalous Hei emission (Szeifert 
et al. I1996f) . NLTE modeling of the complete dataset is 
at present underway to test this assertion; however we 
strongly urge continued monitoring of this star, given the 
motion of W243 across the HD limit into the dynami- 
cally unstable 'yellow void', the rarity of such an event 
and the remarkably high mass rate inferred from current 
observations. In particular, photometric observations to 
determine the presence, or otherwise, of the characteris- 
tic >1 magnitude variability of bona fide LBVs would be 
invaluable. 
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